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| Abstract. Using the narrow clustering of the geometrically corrected gamma-ray 

■ energies released by gamma-ray bursts, we investigate the possibility of using these 
sources as standard candles to probe cosmological parameters such as the matter 
density f2 m and the cosmological constant energy density By simulating different 
samples of gamma-ray bursts, we find that f2 m can be determined with accuracy 
with data from 300 sources. We also show that, if f2 = 1 is due to a quintessence field, 

■ some of the models proposed in the literature may be discriminated from a Universe 
. with cosmological constant, by a similar-sized sample of gamma-ray bursts. 

o 
in 
o 

O ■ 1- Introduction 

ST 

Recent studies have pointed out that Gamma-Ray Bursts (GRBs) may be used as 
standard cosmological candles. The prompt 7-ray energy release, when neglect is made 
of the conical geometry of the emission, spans nearly three orders of magnitude, and 
the distribution of the opening angles of the emission, as deduced from the timing of 
the achromatic steepening of the afterglow emission, spans an identically wide range of 
values. However, when the apparently isotropic energy release and the conic opening 
of the emission are combined to infer the intrinsic, true energy release, the resulting 
distribution does not widen, as is expected for uncorrelated data, but shrinks to a very 
well determined value (Frail et al. 2001; Panaitescu & Kumar 2001; Bloom, Frail, & 
Kulkarni 2003), with a remarkably small (one-sided) scattering, corresponding to about 
a factor of 2 in total energy. Similar studies in the X-ray band (Piran et al. 2001; 
Berger, Kulkarni, Sz Frail 2003) have reproduced the same results. 

It is thus very tempting to study to what extent this property of GRBs makes them 
suitable cosmological standard candles. After an early investigation made by Cohen & 
Piran (1997), Schaefer (2003) proposed using two well known correlations of the GRBs 
luminosity (with variability, and with time delay) to the same end, while Dai, Liang, 
& Xu (2004) and Ghirlanda et al. (2004) exploited the recently reported relationship 
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between the beaming-corrected 7-ray energy and the local observer peak energy in GRBs 
(Ghirlanda, Ghisellini, & Lazzati 2004). We instead neglect these three relationships 
and concentrate on the very narrow spread of the true, geometrically corrected energy 
release as a distance indicator, recalling however that its determination for any given 
bursts requires substantially more information than the methods presented by Schaefer, 
and the other authors mentioned above. 

As for the possible variation of ambient density from burst to burst, which may 
widen the distribution of bursts energies, Frail et al. (2001) remarked that this spread is 
already contained in their data sample, and yet the distribution of energy releases is still 
very narrow. If we were somehow able to measure the distribution ambient densities, 
and subtract these from the sample, the distribution of energy releases should narrow 
even more, not widen: in fact, since we obviously expect the two distributions to be 
uncorrelated, we also expect the one resulting from their combination to be wider than 
the intrinsic distribution of energy releases. 

There are at least two respects in which GRBs are better than type la SuperNovae 
(SNIa) as cosmological candles, one in which they are worse, and one in which they 
are probably even. On the one hand, GRBs are easy to find and locate: even 1980s' 
technology allowed BATSE to locate ~1 GRB per day, despite an incompleteness of 
about 1/3, making the build-up of a 300-object database a one-year enterprise, with 
old technology. The launch of the Swift satellite, which took place on 20 November 2004, 
is expected to detect GRBs at about the same rate as BATSE, but with a nearly perfect 
capacity for identifying their redshifts simultaneously with the afterglow observations ||. 
Second, GRBs have been detected out to very high redshifts: even the current sample 
of about 40 objects (Greiner 2004) contains several events with z > 3, with an absolute 
maximum of z = 4.5 for GRB 000131. This should be contrasted with the difficulty of 
locating SNe with z > 1, and the absolute lack of any SN with z > 2. The currently 
observed distribution of GRBs redshifts contains instead 21 events with z > 1 out of a 
total of 39 (see Figure QJ. 

On the other hand, the distribution of luminosities of SNIa is narrower than 
the distribution of GRBs energy releases, corresponding to a magnitude dispersion 
om = 0.18 rather than om = 0.75. However, the two break even (probably) in 
terms of our understanding of the underlying physical reasons for the uniformity of 
the distributions, which is wanting in both cases. 

Thus GRBs may provide a complementary standard candle, out to distances which 
cannot be probed by SNIa, their major limitation being the larger intrinsic scatter of the 
energy release, as compared to the small scatter in peak luminosities of SNIa. It is thus 
important to assess whether this larger scatter still allows GRBs to be used as standard 
candles. To this end, and as a first aim of the paper, we carry out numerical simulations 
of random samples of GRBs, whose energy releases are distributed as found out by Frail 
et al. (2001), to see to what extent global cosmological parameters can be identified by 

|| http: / / swift.gsfc.nasa.gov/docs/swift /proposals/appendix _f.html 
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an arbitrarily large (but within reason) sample of hypothetical observations. 

As a second aim of the paper, we study, also by means of simulations, whether 
the larger redshift range spanned by GRBs, when compared with SNIa, allows us to 
identify specific models for quintessence. If the non-matter component of the overall 
energy density in the Universe were indeed a constant, at z ~ 1 the increase in the 
matter content would dwarf it, and there would be no difference, at larger redshifts, 
between a model with cosmological constant, and one without (we call this the null 
hypothesis) . However, if the cosmological constant is not constant at all, but is provided 
by the new heuristic field called quintessence, one may hope that at least some models 
display evolution of the cosmological distances (luminosity, fluence, angular, and so on), 
which differ substantially from those of the model with cosmological constant. It is thus 
our second aim to study universes with different, simple models for quintessence, to see 
whether GRBs observations may be able to discriminate between them. In other words, 
we study whether GRBs can reject the hypothesis of a constant A. 

We stress that this paper is not aimed at displaying the potential for cosmological 
investigation by any coming satellite, but instead at determining whether the size of 
a realistically obtainable set of data (perhaps to be obtained by means of a dedicated 
satellite) is useful for cosmological studies. We assume that we know, for every burst in 
our sample, the redshift, and the opening-corrected apparent fluence (i.e., the apparent 
luminosity integrated over the burst duration), and that there is no evolution with 
redshift of the bursts intrinsic energy release. We remark that it is not necessary to 
have a complete and homogeneous sample of objects to carry out this exercise, and that 
the precise value of the bursts average energy release is not necessary, because as usual 
in cosmological tests, we are fitting the dependence of the luminosity distance upon 
redshift and cosmological parameters, not its absolute normalization. 

The plan of the paper is as follows. In Section 2, we display a simulation with the 
simple aim of showing the power of a set of 300 GRBs distributed out to large redshifts, 
in rejecting or accepting the presence of a cosmological constant term in the Universe 
density distribution. A test like this would also be useful in practice, since it would 
be completely independent of observations of fluctuations in the Cosmic Microwave 
Background Radiation (CMBR). Then, in Section 3, we assume a = 0, Q m = 1 
cosmology, and test the ability of similar-size sets of GRBs to determine Q m = 1. In 
Section 4, we assume instead Q\ ^ 0, Q m + Q\ — 1, and test the ability of the same 
samples of GRBs to identify the correct values of Vl m . In Section 5, we abandon the 
hypothesis that A is a constant, and turn to different quintessence models, showing that 
at least one of the important ones (Gasperini, Piazza, & Veneziano 2002) can be easily 
discriminated from the others, and from the null hypothesis. In Section 6, we summarize 
and conclude. 



Can gamma-ray bursts constrain quintessence? 



4 



2. A simple test 

First, in order to show what we are aiming at, we performed a Kolmogorov-Smirnov 
(KS) test on two data sets made of 300 GRBs simulated in two different cosmological 
models, one with fl m = 1 and Qa = and the other with fl m = 0.3 and fl\ = 0.7, but 
both with a Hubble constant H = 65 km s _1 Mpc -1 (as it will be assumed throughout 
the paper). We preferred the KS test to others since it is applicable to any kind of 
continuous distribution that is a function of a single independent variable, which is the 
case we are dealing with. The \ 2 test, for comparison, is more suited to point out 
differences between binned distributions. For the KS test each list of data points, after 
ordering, is converted to a cumulative distribution function giving the fraction of data 
points to the left of a given value for the variable. Then the maximum value D of the 
absolute difference between these two cumulative distribution functions is adopted as 
the test statistic, and the probability Qks of finding values greater than the observed 
D gives the significance level for the null hypothesis that the data sets are drawn from 
the same distribution (Von Mises 1964). 

We assume that GRBs are indeed standard candles with true 7-ray energy released, 
E 1 , following a Gaussian distribution in its logarithm with mean /j, — 51.1 (if E 1 is 
expressed in erg units, Bloom, Frail, & Kulkarni 2003) and o = 0.3 (corresponding to 
a multiplicative factor of 2), and that they are distributed in the Universe according to 
the model of star formation rate Rsfi{z) reported in Porciani & Madau (2001), which 
matches the log A — log P relation (GRB number counts vs. peak photon flux) obtained 
with BATSE data. Applying the KS test to the redshift distributions, we found that 
the probability that the two data sets are the same is Q KS = 0.031, a "no man's land" 
value for this test. On the other hand, the application of the KS test to the parameter 
logd 2 L (z), where di(z) is the luminosity distance, resulted in a significant probability 
Qks ~ 10 -14 , which tells us that it is possible to discriminate between the two different 
cosmological models if a set of 300 GRB luminosity distances is known, without any 
reference to CMBR data. 

3. Simulations in a A = cosmology 

We consider now a A = cosmology, in which the only contribution to the density 
parameter is given by fl m . We assume for GRBs the same energy distribution as for 
the KS test. However, the assumed mean value is not relevant for our investigation, 
since it is the dispersion value that constrains the cosmological density parameter. The 
dispersion of the 7-ray energy released in GRBs may be pinned down in the future by 
a local sample of sources, such as the recently discovered GRB 030329 and 031203, at 
z = 0.1685 (Price et al. 2003) and z = 0.1055 (Malesani et al. 2004) respectively (see 
Section 5 for a more detailed discussion on this point). 

The standard candle energy is related to the fluence of the burst / 7 = E 7 (l + 
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z)/(47rd|(z)) via the luminosity distance d^z), whose expression for A = is: 



dL(Z) = Ho W m « 

Since the /c-correction is independent of any cosmological parameter, we take no account 
of it. A discussion about its effects on the distribution of GRBs energy releases is made 
in Bloom, Frail, & Sari (2001). In order to have a linear propagation of errors throughout 
our simulations, we choose to construct with GRBs a Hubble diagram logd^ — z, since 
the distribution of the parameter \ogd\ is the same of that of log.E 7 , and therefore it is 
Gaussian. 

The number of GRBs per redshift unit is given by the expression: 
dN GRB _ n(z)(dV / dz) 
dz ~ $Q max n(z)(dV / dz)dz 
where n(z) is the redshift distribution function, extending to the maximum redshift for 
GRB explosions z ma x, and dV/dz is the comoving volume element, which for A = is: 

dV = _c ^d 2 L (z) 

dz H (i + z f [Q m (i + z f + (i _ Q m )(i + 2)2]l/2 

As for the KS test, we assume that the redshifts of GRBs are distributed according to 
the model of star formation rate Rsfi( z ) reported in Porciani & Madau (2001): 

exp(3.4s) 

lZjOC exp(3.8z)+45 1 j 

This function increases rapidly between z = and 1, peaks between z = 1 and 2, and 
gently declines at higher redshifts. We fix z max = 5. 

In order to study the ability of GRBs in probing the cosmological parameters as 
a function of their number, we have simulated different samples with Ngrb — 10, 30, 
100, 300 and 1000. Moreover, in order to be free from statistical fluctuations, we have 
performed 10 2 , 10 3 , and 10 4 realizations of each of these samples. 

Now the simulation of a GRB consists of the random sampling of both the redshift z 
and the true 7-ray energy released E 7 , according to the respective adopted distributions. 
Given a cosmological model, from these coupled values we obtain the corresponding 
value for the parameter \ogd 2 L , which we plot on the Hubble diagram as a function of 
z. At this point we perform a x 2 minimization of the simulated data to see with what 
accuracy the fit reproduces the input cosmology. The measurement error on \ogd\ is 
assumed to be a = 0.3. 

In Tabled the mean results of our repeated fits are reported for an input cosmology 
with Q m = 1. First, from this Table it is evident that the mean values obtained 
from the fits are independent of the number of sample realizations, i.e., the intrinsic 
fluctuations corresponding to different samples of GRBs in the same cosmological model, 
are small. Moreover, the Table shows how the accuracy of GRBs in constraining the 
matter density fraction Q m increases with their number Ngrb- The given cosmology is 
readily reproduced by the best fit value for any N GRB , while its dispersion is reduced 
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from ~30% for a sample with N GRB = 10 to ~3% for N GRB = 1000. We have also 
carried out simulations for different values of Q m , but still = 0, obtaining every time 
very similar results. 

It is worth noting that it will be very difficult and time consuming to determine E 1 
of 300 GRBs to the accuracy required. Even then the resulting 7% error on Q m is larger 
than the ~ 1% errors today from WMAP and, eventually, SNAP. Still, an independent 
measurement of a parameter of such paramount importance need not be useless, even 
if late in coming. 

4. Simulations in a A-dominated cosmology 

We move now to a A-dominated cosmology, in which the contributions to the density 
parameter are given by the mass density, fl m , and by the cosmological constant 
energy density, Q\. In the light of the recent observations of the cosmic microwave 
background anisotropy (Bennett et al. 2003), we restrict our simulations to a flat 
Universe Q m + ^a = 1. In this case the expression for the luminosity distance has the 
integral form: 



An analytical fit to this expression, with a relative error of less than 0.4% for 0.2 < 
< 1, is presented in Pen (1999). In order to reduce the run time of our simulations, 

we have exploited this fit to the luminosity distance. 

In a A-dominated cosmology, the number of GRBs per redshift unit is still given 

by equation 0, but in this case the expression for the comoving volume element is: 



For the GRB redshift distribution we adopt the same function as for the A = 
cosmology, with the same value of z max . To take into account the difference in luminosity 
density between an Einstein-de Sitter and a A flat Universe, we applied to n(z) the 
correction factor [f2 m (l + z) 3 + VLt^ 1 ! 2 /{I + z) 3 ^ 2 (see the Appendix of Porciani & Madau 
2001 for details). 

In order to study the ability of GRBs in probing the cosmological parameters in a 
A-dominated Universe, we have simulated 10 2 realizations of GRB samples with Ngrb 
= 10, 30, 100, 300 and 1000. The x 2 minimization of the resulting Hubble diagrams has 
been performed considering \ogd 2 L depending only on the fit parameter Q m , i.e., using 
the relation — 1 — ^m- Table El reports the general results of our repeated fits for a 
flat cosmology with input values Q m = 0.3 and = 0.7 (which are those adopted in 
Frail et al. 2001). 

As in the A = case, the accuracy of GRBs in constraining the two contributions 
Q m and ^a to the density parameter increases with their number Ngrb, reducing the 
dispersion about the best value for the fit parameter Q m from ~40% for a sample with 
N GRB = 10 to ~4% for Ngrb = 1000. 




(5) 



dV c ^d 2 L (z) 



(6) 



dz ^o(l + ^) 2 [fi m (l + z)3 + fi A ] 



1/2 
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Figure 1. Hubble diagram \ogd 2 L — z with data simulated for a sample of 300 GRBs 
in a flat Universe with density parameters Q m — 0.3 and Q\ = 0.7. The solid curve 
shows the function logd^(z) in the assumed cosmology, while the dashed curves give 
the dispersion about the best fit parameter (the upper curve corresponds to lower Q m ). 

Focussing on the samples in a A-dominated cosmology with N GRB = 300, a data set 
which can be realistically obtained in future space missions, Figure Q shows one of the 
Hubble diagrams log d\ — z obtained with the simulations. The solid curve shows the 
function log d\(z) in the assumed cosmology, while the dashed curves give the dispersion 
about the best fit parameters. The ability of a sample of 300 GRBs in constraining the 
actual cosmology is evident. The statistical fluctuations of the Nqrb = 300 sample fit 
are outlined in the histogram of Figure which shows the distribution of the best fit 
values of the matter density fraction Q m for 10 3 sample realizations. The distribution 
peaks at Q m = 0.3, has a dispersion Sa m = 0.0228, and a kurtosis ku m = 3.0993, to be 
compared with the value of a Gaussian distribution, i.e., 3. 

We have also investigated the effects of changing the assumptions of our simulations 
on the probing ability of GRB samples to determine the geometry of the Universe. 
First, we have considered a GRB redshift distribution ruled by the simple function 
n(z) oc (1 + z) 3 instead of equation (jH). The result is a slight decrease of the dispersion 
about the best fit parameter Q m at all values of Ngrb- This is due to the larger number 
of GRBs sampled at high redshift values by this alternative distibution, which increases 
monotonically with z. At high redshifts the distinction between different cosmologies 
becomes more evident (see curves in Figure^), thus more GRBs at large z imply better 
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Figure 2. Histogram with the distribution of the best fit values of the matter density 
fi m for 10 3 realizations of a sample of 300 GRBs in a flat Universe with density 
parameters fl m = 0.3 and Oa = 0.7. The distribution has a mean < fi m >= 0.3001, a 
median Q m (med) — 0.3002, a dispersion Sn m = 0.0228, and a kurtosis ka m = 3.0993. 

constraints on the cosmological parameters. 

Then we have studied the effect of varying the dispersion about the standard candle 
energy. We assumed a = 0.6 (doubling the dispersion to a multiplicative factor of 4) 
about the logarithmic mean value fi — 51.1 reported in Bloom, Frail, & Kulkarni (2003). 
We find that the resulting effect is of course a worse accuracy in the reproduction of the 
input cosmology, the dispersion about the best fit parameter Vt m increasing by a factor 
~ 2 at all values of Ngrb- In particular, it is ~15% for a sample with Ngrb = 300. 

Moreover, we point out that the variation of the standard candle energy E 1 has no 
effect on the ability of GRB samples in putting constraints on cosmological parameters, 
since the mean value of the Gaussian distribution of log E 1 gives only a normalization 
constant to our simulations, but is not instructive for their scatter. 

Finally, we must remark that the analyses of both Frail et al. (2001) and Bloom, 
Frail, & Kulkarni (2003) assume of course a particular set of cosmological parameters 
(Vt m = 0.3, = 0.7, and H = 65 km s -1 Mpc -1 ) to derive the standard 7-ray energy 
of GRBs. To avoid a circular logic and the limitations in cosmographic applications 
pointed out in Bloom, Frail, & Kulkarni (2003), we should assume a candle calibration 
with a local sample of sources, a prospect which can now be considered possible in the 
light of the discovery of the nearby GRBs 030329 and 031203 (see Section 5 for a more 
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detailed discussion on this point). 

5. Simulations in a quintessence cosmology 

Now we abandon the cosmological constant Universe and we consider some of the most 
popular quintessence models. In particular, we choose as tracker potential classes the 
inverse power-law Ratra-Peebles potential (hereafter RP; Ratra & Peebles 1988), defined 
as: 



Following Caresia, Matarrese, & Moscardini (2004), we considered such potentials within 
the framework of "extended quintessence" models. These are characterized by a coupling 
between gravity and quintessence ruled by a parameter £, where £ = means no 
coupling. By fixing a and £, we obtained numerical values for the luminosity distance 
cIl(z) with the standard procedure. First, we numerically solved the Klein-Gordon 
equation for the scalar field 0. From this solution it is possible to get the expansion rate 
H(z) = H h(z), which depends on <p since the quintessence scalar field contributes to 
the total energy density. Then we exploited the usual relation between expansion rate 
and luminosity distance: 



We repeated the same procedure for the dilaton scenario introduced by Gasperini, 
Piazza, & Veneziano (2002; hereafter GPV). In order to distinguish between different 
tracker potentials of the same class, say RP, we will use a couple of indices, the first 
one giving the value of the a parameter, while the second one referring to the adopted 
coupling parameter £, being 1 for £ = and 2 for £ = 0.01 (therefore, model RP01 
corresponds to no quintessence). In Figure E3 we report the function \ogd 2 L (z) found 
for the two most "extreme" among the quintessence models considered, i.e., those in 
which this function differs most strongly from that corresponding to a A-dominated flat 
Universe with density parameters Q m = 0.3 and = 0.7, already shown in Figure ^ 
All the other models give a \og d 2 L (z) curve lying between the no quintessence case and 
the RP22 curve. 

Figure 01 allows us to make an important, technical point. We have so far postponed 
the thorny issue of the calibration of the absolute energy release: in fact, when a 
sample of GRBs is used to derive the distribution of energy releases, it is necessary 
to assume a set of cosmological parameters to compare GRBs at different redshifts, 
making (potentially!) our argument circular. However, it is well known that, for z <C 1, 
all cosmological models coincide, so that a local calibration of the absolute energy release 




(7) 



and the SUGRA potential (Binetruy 1999; Brax & Martin 1999): 




(8) 




(9) 
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Figure 3. The parameter logd^ as a function of redshift z in a flat Universe with 
density parameters f2 m = 0.3 and Qa = 0.7 (solid curve) compared with those obtained 
in the quintessence models RP22 (lower dashed curve) and GPV (upper dotted curve). 



is possible. Originally, given the very large redshifts of the first GRBs, it was not clear 
whether this could be achieved, but with the accumulation of further data this does not 
look like a real concern: there are currently 4 GRBs with z < 0.4, out of a total of 39 
(Figure EJ), making a local calibration with a large sample of GRBs a real possibility. 

Furthermore, the situation is even better if we assume f2 m -\- £l\ — 1: we see from 
Figure El that all models, including GPV which is by far the most discordant one, yield 
essentially the same luminosity distance out to z « 1. There are currently 17 GRBs 
with z < 1 out of a total of 39 (Figure 0]), making the issue of calibration a moot one, 
once a sufficiently large sample is obtained. 

In order to determine the number of GRBs per redshift unit in quintessence 
universes, we can exploit again the Appendix of Porciani & Madau (2001) and find 
that: 

dN GRB dj(z)n(z)/(l + zy/i 

alz f* ma *d 2 L (z)n(z)/(l + z) 7 / 2 dz 1 ' 

independently of dV/dz, where we have already taken into account the correction factor 
for the difference in luminosity density from the Einstein-de Sitter Universe. The 
adopted GRB redshift distribution n(z) is always given by equation (HJ), and we fix 
again z max = 5. 
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Our purpose is now to investigate whether it is possible to discriminate between 
different quintessence cosmological models via a set of GRBs, considered as standard 
candles. To do this, we performed a series of KS tests on simulated data sets with Nqrb 
= 100, 300 and 1000 in two different A-dominated flat cosmologies both with Q m = 0.3 
and fix = 0.7 at z — 0, but, while one has a truly constant A, the other one has a 
quintessence field which reduces to = 0.7 at z = 0. The model for quintessence is 
chosen from those listed before. We considered both the case in which the same set 
of GRB redshifts is used, and the case in which different patches of the Universe are 
selected, each with its own cosmology, and thus two fully distinct sets of GRB redshifts 
are used for the two different cosmologies. In this second way, we wish to include the 
cosmic variance into our simulations. Our results are presented in Table El for a same 
GRB redshift distribution, and in Table EJ where cosmic variance has been considered. 

From Tables El and |H we see that it is quite difficult to discriminate a quintessence 
cosmological model using a set of up to 1000 GRBs as standard candles. Only the GPV 
model could be significantly discriminated with a set of 1000 GRBs, especially if we 
take into account cosmic variance, in which case we obtained < Qks >— 4.89 ■ 10~ 8 . 
A discrimination may also be possible for the RP22 model, which resulted in a 
< Qks >= 8.45 • 10~ 3 with 1000 GRBs, when not considering cosmic variance. These 
results are consistent with the extreme behaviours of the log d 2 L (z) functions for GPV 
and RP22 models, as already stated before and shown in Figure El 

These findings may be the consequence of the high dispersion, a = 0.3 in its 
logarithm, around the mean value of the geometrically corrected 7-ray energy released 
by GRBs, which we have taken as a cosmological candle. In order to investigate whether 
a better discrimination of different quintessence models could be obtained with a pin 
down of our standard candle, we performed another series of KS tests on simulated data 
sets, in which we adopted a SNIa-like candle with a logarithmic dispersion of a = 0.072, 
corresponding to a conservative magnitude dispersion of o~m = 0.18 mag, intermediate 
between the o~m = 0.21 mag value found by Riess, Press, & Kirshner (1995) and the 
later om = 0.12 mag of the same authors (Riess, Press, & Kirshner 1996) together with 
the very recent results by Wang et al. (2003) in the range o~m = 0.08 — 0.11 mag. 
The conditions of the simulations were the same as for the first KS test series, with 
the additional subdivision of a value z max = 5, typical of GRBs, or 2, more typical of 
SNIa. In the first case, our results are reported in Table for the same standard candle 
redshift distribution, and in Table EH where cosmic variance has been considered, while 
in the second case the results are shown respectively in Tableland Table El 

Analysis of Tables EHH1 shows that with a less dispersed cosmological candle the 
GPV and RP22 quintessence models could be significantly discriminated with a set of 
1000 sources observed up to z max = 5, the first one in this case perhaps even with only 
300 sources (< Qks > = 1-04 • 10~ 3 taking into account cosmic variance). Moreover, 
a hint for discrimination is also given with 1000 candles by the RP21 model, which 
resulted in a < Qks >— 6.92 • 10 -3 when not considering cosmic variance. On the other 
hand, limiting the candle distribution to z max = 2, it seems possible to significantly 
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discriminate all the quintessence models of the RP tracker potential class, together with 
the GPV model, if 1000 sources were observed. The RP21 and RP22 models may 
be discriminated even with a data set of only 300 standard candles. These results 
are consistent with Figure El where it is possible to notice how up to z = 2 the 
logd 2 L (z) function for the RP22 (and also for the not reported RP21) model differs 
more than the GPV one from the no quintessence case. There seems to be no chance for 
a discrimination of the SUGRA class of quintessence models, whose \ogd\(z) functions 
are in fact very similar to the solide curve of Figure El 

6. Conclusions 

We have simulated different samples of GRBs adopting 7-ray energy and redshift 
distributions consistent with recent observational results, in order to investigate their 
ability to probe cosmological parameters such as the density fractions Q m and Q\. Our 
result is that in a A-dominated flat Universe the accuracy in the determination of the 
matter density fl m is ~40% for a sample with Ngrb = 10 and an excellent ~4% for 
N GRB = 1000. 

For comparison, a ~20% accuracy on the determination of fl m has been recently 
claimed by using GRBs as standard candles ruled by the luminosity-variability and 
luminosity-lag time relations (Takahashi et al. 2003). 

Since GRBs are much more readily observed than SNIa, especially at high redshifts 
(notice that in the sample observed so far, which is reported in Figure 01 there are 
already 4 GRBs with z > 3, about 10% of the whole), they should allow us to probe 
cosmological parameters more deeply than these latter sources. Moreover, during the 
last few years of observations the number of GRBs with known redshifts has almost 
reached the same number of high redshift SNIa discovered by the Supernova Cosmology 
Project (Perlmutter et al. 1999). 

Lastly, after showing that the absolute energy release can be calibrated using the 
low redshift GRBs, we have shown that GRBs have the potential to investigate the 
luminosity distance out to large redshifts, and this, in turn, means that at least some 
models for quintessence, among which the important dilaton model of Gasperini, Piazza 
and Veneziano, can be tested and discriminated from competing models. 

On 20 November 2004, the Swift satellite was launched, and the detection of ~200 
GRBs with known redshifts is expected during the 3 years of its sky observations. If 
GRBs are confirmed by these new data to emit a standard amount of energy, then our 
simulations stress how the accuracy in the determination of cosmological parameters 
increases with the number of their known redshifts, making Swift, at least potentially, 
a GRB Cosmology Project. 
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Figure 4. Distribution of all the known redshifts of GRBs as of November 2004. All 
redshifts are within the 0.1—4.5 range, apart from GRB 980425, possibly associated 
with the nearby SN1998bw (z=0.0085, Galama et al. 1998). Data is taken from Greiner 
(2004). 
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Table 1. Mean values of the fitted cosmological density parameter fl m , of its 
error Af2 m and its dispersion Sn m obtained by fitting 10 2 (top), 10 3 (middle) and 
10 4 (bottom) GRB sample realizations with Nqrb distributed according to function 
Rsfi{z) of Porciani & Madau (2001) in an Einstein-de Sitter Universe (fi TO = 1). 



Nqrb 


< > 


< AO m > 




10 


0.9983 


0.2997 


0.3097 


30 


1.0158 


0.1895 


0.1993 


100 


0.9937 


0.0993 


0.1108 


300 


0.9959 


0.0599 


0.0629 


1000 


1.0009 


0.0332 


0.0351 



Nqrb 


<n m > 


< AO m > 




10 


1.0265 


0.3095 


0.3483 


30 


1.0142 


0.1766 


0.1993 


100 


1.0026 


0.0993 


0.1048 


300 


1.0020 


0.0593 


0.0595 


1000 


1.0015 


0.0331 


0.0333 



Ngrb 




< AO m > 




10 


1.0320 


0.3085 


0.3516 


30 


1.0103 


0.1780 


0.1939 


100 


1.0043 


0.0997 


0.1074 


300 


1.0011 


0.0590 


0.0614 


1000 


1.0004 


0.0330 


0.0332 



Table 2. Mean values of the fitted cosmological density parameters TO and Q\, of 
their error Afi and their dispersion Sq obtained by fitting 10 2 GRB sample realizations 
with Nqrb distributed according to function Rsfi(z) of Porciani & Madau (2001) in 
a flat Universe with input values Q m = 0.3 and Oa = 0.7. 



Nqrb 


< n m > 


< > 


< Afl > 


Sn 


10 


0.3195 


0.6805 


0.1004 


0.1307 


30 


0.2973 


0.7027 


0.0763 


0.0700 


100 


0.3002 


0.6998 


0.0363 


0.0351 


300 


0.3023 


0.6977 


0.0219 


0.0222 


1000 


0.3001 


0.6999 


0.0120 


0.0125 
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Table 3. Mean values of the probability Qks and the maximum distance D for 
a KS test on the parameter \ogd 2 L (z) of 100 realizations of a redshift sample made 
of 100 (top), 300 (middle) and 1000 (bottom) GRBs obtained in a A-dominatcd fiat 
cosmology with TO = 0.3 and Oa = 0.7 at z = 0, but with two different logd|(z) 
distributions, one resulting from a cosmology with a truly constant A and the other 
from a quintessence model defined by its name in the first column, which reduces to 
the same cosmology at z = 0. 



Model 


< Q 


KS 


> 


< D > 


RP01 


6.84 


■ 10 


-1 


0.098 


RP11 


6.08 


• 10 


-1 


0.107 


RP12 


5.75 


• 10 


-1 


0.110 


RP21 


4.91 


• 10 


-1 


0.120 


RP22 


4.54 


• 10 


-1 


0.125 


SUGRA11 


7.29 


• 10 


-1 


0.094 


SUGRA12 


7.25 


• 10 


-1 


0.094 


SUGRA21 


6.37 


• 10 


-1 


0.104 


SUGRA22 


6.52 


• 10 


-1 


0.103 


GPV 


2.76 


• 10 


-1 


0.149 




Model 


< Q 


KS 


> 


< D > 


RP01 


6.98 


■ 10 


-1 


0.057 


RP11 


4.37 


• 10 


-1 


0.073 


RP12 


4.22 


• 10 


-1 


0.076 


RP21 


2.14 


• 10 


-1 


0.094 


RP22 


1.41 


• 10 


-1 


0.102 


SUGRA11 


6.38 


• 10 


-1 


0.060 


SUGRA12 


6.12 


• 10 


-1 


0.061 


SUGRA21 


5.05 


• 10 


-1 


0.069 


SUGRA22 


5.26 


• 10 


-1 


0.067 


GPV 


8.32 


• 10 


-2 


0.113 



Model 


< Q 


KS 


> 


< D > 


RP01 


7.16 


10 


-1 


0.031 


RP11 


1.85 


10 


-1 


0.053 


RP12 


1.35 


10 


-1 


0.056 


RP21 


1.98 


10 


-2 


0.076 


RP22 


8.45 


10 


-3 


0.082 


SUGRA11 


4.70 


10 


-1 


0.039 


SUGRA12 


4.56 


10 


-1 


0.040 


SUGRA21 


3.19 


10 


-1 


0.045 


SUGRA22 


2.94 


10 


-1 


0.047 


GPV 


3.28 


10 


-4 


0.103 
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Table 4. Mean values of the probability Qks and the maximum distance D for a KS 
test on the parameter log d\(z) of 100 realizations of two different samples made of 
100 (top), 300 (middle) and 1000 (bottom) GRBs, one obtained in a A-dominated flat 
cosmology with Q m — 0.3 and Q\ = 0.7, and the other in a quintessence model defined 
by its name in the first column, which reduces to the same cosmology at z = 0. 



Model 


< Q 


KS 


> 


< D > 


RP01 


4.64 


10 


-1 


0.126 


RP11 


4.91 


10 


-1 


0.123 


RP12 


4.91 


10 


-1 


0.123 


RP21 


4.60 


10 


-1 


0.128 


RP22 


3.55 


10 


-1 


0.143 


SUGRA11 


4.97 


10 


-1 


0.121 


SUGRA12 


4.92 


10 


-1 


0.121 


SUGRA21 


4.52 


10 


-1 


0.127 


SUGRA22 


5.42 


10 


-1 


0.117 


GPV 


1.41 


10 


-1 


0.190 



Model 


< Q 


KS 


> 


< D > 


RP01 


5.01 


10 


-1 


0.070 


RP11 


3.45 


10 


-1 


0.082 


RP12 


4.16 


10 


-1 


0.076 


RP21 


2.30 


10 


-1 


0.099 


RP22 


1.80 


10 


-1 


0.106 


SUGRA11 


5.16 


10 


-1 


0.069 


SUGRA12 


5.06 


10 


-1 


0.070 


SUGRA21 


4.18 


10 


-1 


0.077 


SUGRA22 


4.50 


10 


-1 


0.074 


GPV 


7.36 


10 


-3 


0.173 




Model 


< Q 


KS 


> 


< D > 


RP01 


5.31 


10 


-1 


0.037 


RP11 


2.32 


10 


-1 


0.053 


RP12 


1.49 


10 


-1 


0.062 


RP21 


5.43 


10 


-2 


0.077 


RP22 


1.81 


10 


-2 


0.086 


SUGRA11 


3.74 


10 


-1 


0.045 


SUGRA12 


4.09 


10 


-1 


0.043 


SUGRA21 


3.54 


10 


-1 


0.047 


SUGRA22 


2.66 


10 


-1 


0.051 


GPV 


4.89 


10 


-8 


0.158 
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Table 5. The same as Table 01 but obtained by adopting SNIa-like standard candles 
distributed with a logarithmic dispersion a — 0.072 up to z max — 5. 



Model 


< Q 


KS 


> 


< D > 


RP01 


9.66 


10 


_ I 


0.062 


RP11 


7.83 


10 


-1 


0.090 


RP12 


7.20 


10 


-1 


0.096 


RP21 


5.45 


10 


-1 


0.112 


RP22 


4.27 


10 


-1 


0.125 


SUGRA11 


9.21 


10 


-1 


0.070 


SUGRA12 


8.99 


10 


-1 


0.075 


SUGRA21 


8.52 


10 


-1 


0.081 


SUGRA22 


8.21 


10 


-1 


0.086 


GPV 


1.93 


10 


-1 


0.157 



Model 


< Q 


KS 


> 


< D > 


RP01 


9.68 


10 


-1 


0.037 


RP11 


5.11 


10 


-1 


0.067 


RP12 


4.45 


10 


-1 


0.071 


RP21 


1.65 


10 


-1 


0.092 


RP22 


1.09 


10 


-1 


0.099 


SUGRA11 


8.19 


10 


-1 


0.050 


SUGRA12 


8.18 


10 


-1 


0.050 


SUGRA21 


6.59 


10 


-1 


0.059 


SUGRA22 


6.39 


10 


-1 


0.060 


GPV 


2.03 


10 


-2 


0.129 



Model 


< Q 


KS 


> 


< D > 


RP01 


9.68 


■ 10 


-1 


0.020 


RP11 


1.39 


■ 10 


-1 


0.052 


RP12 


9.51 


• 10 


-2 


0.056 


RP21 


6.92 


• 10 


-3 


0.076 


RP22 


1.54 


■ 10 


-3 


0.087 


SUGRA11 


5.49 


■ 10 


-1 


0.036 


SUGRA12 


5.17 


■ 10 


-1 


0.037 


SUGRA21 


2.95 


■ 10 


-1 


0.044 


SUGRA22 


2.45 


■ 10 


-1 


0.046 


GPV 


5.57 


• 10 


-6 


0.118 
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Table 6. The same as Table 0] but obtained by adopting SNIa-like standard candles 
distributed with a logarithmic dispersion a — 0.072 up to z max — 5. 



Model 


< Q 


KS 


> 


< D > 


RP01 


5.25 


10 


_ i 


0.118 


RP11 


4.81 


10 


-l 


0.122 


RP12 


4.55 


10 


-l 


0.127 


RP21 


4.06 


10 


-l 


0.138 


RP22 


3.36 


10 


-l 


0.145 


SUGRA11 


4.59 


10 


-i 


0.126 


SUGRA12 


4.81 


10 


-l 


0.121 


SUGRA21 


4.66 


10 


-l 


0.125 


SUGRA22 


4.68 


10 


-l 


0.128 


GPV 


6.28 


10 


-2 


0.217 



Model 


< Q 


KS 


> 


< D > 


RP01 


5.07 


10 


-1 


0.069 


RP11 


3.56 


10 


-1 


0.083 


RP12 


3.21 


10 


-1 


0.088 


RP21 


1.97 


10 


-1 


0.101 


RP22 


1.09 


10 


-1 


0.118 


SUGRA11 


4.52 


10 


-1 


0.075 


SUGRA12 


3.92 


10 


-1 


0.079 


SUGRA21 


4.14 


10 


-1 


0.078 


SUGRA22 


3.38 


10 


-1 


0.083 


GPV 


1.04 


10 


-3 


0.191 



Model 


< Q 


KS 


> 


< D > 


RP01 


5.22 


10 


-1 


0.038 


RP11 


1.57 


10 


-1 


0.059 


RP12 


1.27 


10 


-1 


0.063 


RP21 


1.89 


10 


-2 


0.082 


RP22 


1.58 


10 


-2 


0.089 


SUGRA11 


3.58 


10 


-1 


0.045 


SUGRA12 


3.62 


10 


-1 


0.046 


SUGRA21 


2.27 


10 


-1 


0.053 


SUGRA22 


2.43 


10 


-1 


0.051 


GPV 


6.07- 


10" 


-10 


0.177 
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Table 7. The same as Table 01 but obtained by adopting SNIa-like standard candles 
distributed with a logarithmic dispersion a — 0.072 up to z max — 2. 



Model 


< Q 


KS 


> 


< D > 


RP01 


8.78 


10 


_ I 


0.077 


RP11 


3.75 


10 


-1 


0.131 


RP12 


2.99 


10 


-1 


0.141 


RP21 


8.10 


10 


-2 


0.188 


RP22 


6.73 


10 


-2 


0.192 


SUGRA11 


7.14 


10 


-1 


0.096 


SUGRA12 


6.61 


10 


-1 


0.101 


SUGRA21 


5.55 


10 


-1 


0.113 


SUGRA22 


5.12 


10 


-1 


0.116 


GPV 


3.33 


10 


-1 


0.140 



Model 


< Q 


KS 


> 


< D > 


RP01 


8.89 


10 


-1 


0.044 


RP11 


7.51 


10 


-2 


0.108 


RP12 


5.24 


10 


-2 


0.116 


RP21 


2.32 


10 


-3 


0.157 


RP22 


5.57 


10 


-4 


0.176 


SUGRA11 


4.56 


10 


-1 


0.071 


SUGRA12 


4.33 


10 


-1 


0.072 


SUGRA21 


2.19 


10 


-1 


0.089 


SUGRA22 


1.80 


10 


-1 


0.092 


GPV 


4.93 


10 


-2 


0.116 



Model 


< Q 


KS 


> 


< D > 


RP01 


8.77 


■ 10 


-1 


0.025 


RP11 


1.03 


■ 10 


-3 


0.094 


RP12 


2.90 


• 10 


-4 


0.100 


RP21 


5.57 


■ 10 


-8 


0.143 


RP22 


1.53 


■ 10 


-9 


0.158 


SUGRA11 


1.10 


■ 10 


-1 


0.057 


SUGRA12 


1.04 


■ 10 


-1 


0.057 


SUGRA21 


1.85 


■ 10 


-2 


0.072 


SUGRA22 


1.22 


■ 10 


-2 


0.076 


GPV 


1.25 


• 10 


-4 


0.104 
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Table 8. The same as Table 0] but obtained by adopting SNIa-like standard candles 
distributed with a logarithmic dispersion a — 0.072 up to z max — 2. 



Model 


< Q 


KS 


> 


< D > 


RP01 


4.94 


10 


_ I 


0.122 


RP11 


2.89 


10 


-1 


0.153 


RP12 


3.09 


10 


-1 


0.150 


RP21 


1.50 


10 


-1 


0.189 


RP22 


1.01 


10 


-1 


0.199 


SUGRA11 


4.53 


10 


-1 


0.128 


SUGRA12 


3.92 


10 


-1 


0.135 


SUGRA21 


3.82 


10 


-1 


0.139 


SUGRA22 


3.52 


10 


-1 


0.142 


GPV 


2.42 


10 


-1 


0.166 



Model 


< Q 


KS 


> 


< D > 


RP01 


4.81 


10 


-1 


0.070 


RP11 


1.23 


10 


-1 


0.116 


RP12 


8.68 


10 


-2 


0.117 


RP21 


8.56 


10 


-3 


0.162 


RP22 


2.67 


10 


-3 


0.179 


SUGRA11 


3.27 


10 


-1 


0.084 


SUGRA12 


3.41 


10 


-1 


0.084 


SUGRA21 


2.30 


10 


-1 


0.097 


SUGRA22 


2.19 


10 


-1 


0.098 


GPV 


3.89 


10 


-2 


0.138 



Model 


< Q 


KS 


> 


< D > 


RP01 


5.38 


■ 10 


-1 


0.037 


RP11 


6.83 


■ 10 


-3 


0.097 


RP12 


1.11 


• 10 


-3 


0.103 


RP21 


1.32 


■ 10 


-7 


0.149 


RP22 


1.63 


■ 10 


-8 


0.164 


SUGRA11 


1.44 


■ 10 


-1 


0.059 


SUGRA12 


1.32 


■ 10 


-1 


0.061 


SUGRA21 


3.56 


■ 10 


-2 


0.077 


SUGRA22 


3.49 


■ 10 


-2 


0.079 


GPV 


5.27 


• 10 


-5 


0.124 



